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7— ray bursts (GRBs) are likely sources of ultra-high energy, > 10 eV, protons and high energy, 
> 1 TeV, neutrinos. Large volume detectors of ultra high energy cosmic rays (UHECRs) and high 
energy neutrinos, which are already operating and are being expanded, may allow to test in the 
coming few years the predictions of the GRB model for high energy proton and neutrino production. 
Detection of the predicted signals will allow to identify the sources of UHECRs and will provide a 
unique probe, which may allow to resolve open questions related to the underlying physics of GRB 
1 models. Moreover, detection of GRB neutrinos will allow to test for neutrino properties (e.g., flavor 

oscillations for which t 's would be a unique signature, and coupling to gravity) with an accuracy 
many orders of magnitude better than is currently possible. 
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The cosmic-ray spectrum extends to energies ~ 10 eV, and is likely dominated beyond 
~ 10 19 eV by extra-Galactic sources of protons [1] (see, however, [2D). The origin of the 
highest energy, > 10 19 eV, cosmic rays (UHECRs) is a mystery. As explained in §12.1[ the 
stringent constraints, which are imposed on the properties of possible UHECR sources 
by the high energies observed, rule out almost all source candidates, and suggest that 
y-ray bursts (GRBs) and active galactic nuclei (AGN) are the most plausible sources. 
In § 12.21 we show that the energy loss of UHECR protons due to interaction with 
microwave background photons, which limits the proton propagation distance, disfavors 
AGNs as UHECR sources. It is also shown that both the energy production rate and 
the spectrum of UHECRs are consistent with those expected in the GRB model for 
UHECR production. Unique predictions of the GRB model, which may be tested with 
large UHECR detectors, are discussed in § 12.31 (for a pedagogical review of UHECR 
production in GRBs see [3Q). 

High energy neutrino production is discussed in § O It is first shown in § 13.11 that 
cosmic ray observations set an upper limit of Ey<& v < Ey^f B = 5 x 10~ 8 GeV/cm 2 ssr 
on the diffuse extra-Galactic high energy neutrino intensity produced by sources which, 
like GRBs and AGN jets, are optically thin for high-energy nucleons to py and pp(n) 
interactions. The implications of the upper bound to the detector size required for 
detecting extra-Galactic high energy neutrinos are briefly discussed. In §13.21 we show 
that production of 100 TeV neutrinos in the region where GRB y-rays are produced 
is a genericprediction of the GRB fireball model (for detailed reviews of the fireball 
model see [4]; for a pedagogical review of high energy neutrino production in GRBs 
see [H]). 100 TeV neutrino production is a direct consequence of the assumptions that 



energy is carried from the underlying engine, most likely a (few) solar mass black hole, 
as kinetic energy of protons and that y-rays are produced by synchrotron emission of 
shock accelerated electrons. The detection of the predicted neutrino signal will therefore 
provide strong support for the validity of underlying model assumptions, which is 
difficult to obtain using photon observations (due to the high optical depth in the vicinity 
of the GRB "engine"). The predicted neutrino intensity, « 0.2<&^ B , implies a detection 
of ~ 20 muon induced neutrino events per yr in a km-scale neutrino detector. Since 
these events should be correlated in time and direction with GRB y-rays, the search for 
GRB neutrinos is essentially background free. Neutrinos may be produced also in other 
stages of fireball evolution, at energies different than 100 TeV. The production of these 
neutrinos is dependent on additional model assumptions. As an example, we discuss in 
§13.3l the production of TeV neutrinos expected in the "collapsar" scenario, where GRB 
progenitors are associated with the collapse of massive stars. 

The discussion of GRB neutrino emission demonstrates that in addition to identifying 
the sources of UHECRs, high-energy neutrino telescopes can also provide a unique 
probe of the physics of these sources. Moreover, detection of high energy neutrinos 
from GRBs may also provide information on fundamental neutrino properties [5]. High 
energy neutrinos are expected to be produced in astrophysical sources by the decay of 
charged pions, which lead to the production of neutrinos with flavor ratio <3> V( , : ^v M : 
4> Vt = 1:2:0 (here 4> V/ stands for the combined flux of V/ and V/). Neutrino oscillations 
then lead to an observed flux ratio on Earth of 4> V( , : <J> Vm : ^>v T = 1:1:1 [6]. Up- 
going t's, rather than /i's, would be a distinctive signature of such oscillations, and 
it has been pointed out |a 0] that searching for deviations from the standard flavor 
ratio 1:1:1 may enable one to probe new physics. Detection of neutrinos from GRBs 
could, moreover, be used to test the simultaneity of neutrino and photon arrival to 
an accuracy of ~ 1 s, checking the assumption of special relativity that photons and 
neutrinos have the same limiting speed. These observations would also test the weak 
equivalence principle, according to which photons and neutrinos should suffer the same 
time delay as they pass through a gravitational potential. With 1 s accuracy, a burst at 
1 Gpc would reveal a fractional difference in limiting speed of 10~ 17 , and a fractional 
difference in gravitational time delay of order 10~ 6 (considering the Galactic potential 
alone). Previous applications of these ideas to supernova 1987 A (see [|80 for review), 
yielded much weaker upper limits: of order 10~ 8 and 10~ 2 respectively. 

In order to illustrate the possibilities that may be opened by high energy neutrino 
observations, we show in § 13.41 that electromagnetic energy losses of tt's and /i's 
modify the flavor ratio (measured at Earth) of neutrinos produced by it decay, 4> Vc : 
4> Vfl : 4> Vt , from 1 : 1 : 1 at low energy to 1 : 1.8 : 1.8 at high energy. For GRBs the 
transition is expected at ~ 100 TeV, and may be detected by km-scale v telescopes. 
While the detection of x neutrinos would test flavor oscillations, a measurement of the 
flavor transition energy and energy-width will provide a unique probe of the physical 
conditions in the y-ray production region. 



2. ULTRA HIGH ENERGY COSMIC RAYS 



2.1. The acceleration challenge 

Most models of particle acceleration in astrophysical sources involve the acceleration 
of charged particles by the electro-motive force produced through the motion of a 
magnetized plasma. General phenomenological considerations imply that, independent 
of the details of the acceleration mechanism, the minimum energy output from a source 
capable of accelerating in this manner a proton to energy E p is |9|] 

^fd) 2 - 1045 - 5 ?^) 2 ^ <•> 

Here, u = /3c is the characteristic plasma velocity, and T = (1 — /3 2 )~ 1//2 . Only two 
types of sources are known to satisfy this requirement. The brightest steady sources 
are active galactic nuclei (AGN). For them T is typically between 3 and 10, implying 
L > 10 47 erg/ s, which may be satisfied by the brightest AGN fl 1 ON . The brightest transient 
sources are GRBs. For these sources T ~ 10 25 implying L > 10 50 5 erg/s, which is 
generally satisfied since the typical observed MeV-photon luminosity of these sources is 
L r ~ 10 52 erg/s. 

The constraint of eq. © is essentially obtained by requiring the acceleration time to be 
shorter than the proton confinement time. A second constraint is imposed by requiring 
the proton acceleration time to be smaller than its energy loss time. For a relativistic 
wind, where acceleration takes place in internal shocks arising from variability of the 
underlying source driving the wind on time scale At, the proton acceleration time is 



smaller than its synchrotron energy loss time provided 111 111 
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This constraint is remarkably similar to that inferred from y-ray observations on the 
relativistic dissipative winds which are assumed, within the context of the fireball model 
(see [0] for reviews), to produce GRBs: T > 300 is implied by the y-ray spectrum and 
by the short variability time, At ~ 10ms, through the requirement to avoid high pair- 
production optical depth. This, combined with the constraint of eq. was one of the 
two main arguments suggested in 111 ill in support of an association between UHECR and 
GRB sources. 



2.2. The GZK effect 

The energy loss of UHECR protons due to interaction with microwave background 
photons (the "GZK effect") limits the proton propagation distance to < 100 Mpc for 
E p > 10 20 eV OJ]. This disfavors AGNs as UHECR sources, since no AGN powerful 
enough to satisfy the constraint of eq. ([]]) is known to lie within 100 Mpc from Earth. 



Since most GRBs are known to reside at redshift z > 1, it may appear that they too are 
disfavored by the short, 100 Mpc, propagation distance of the protons. This is, however, 
not the case. While GRB y-rays arrive as a burst, slight deflections of protons by an 
inter-galactic magnetic field would cause large delays in their arrival time (compared to 
the photon arrival time). The arrival time delay is accompanied by a comparable arrival 
time spread flll|1 . 
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where D is the source distance, B and X are the field strength and correlation length, 
and a limit 5A 1//2 < 10 8 G Mpc 1 / 2 is set by Faraday rotation measurements lfl3n. The 
number of GRBs contributing to the flux above energy E p at any given time is lllll. [l4ll 
AW> E p ) « (47t/5)R GRB D 3 c {Ep)T[E p ,D c (Ep)], where R GRB is the local (z = 0) GRB 
rate density and D c (E p ) is the proton propagation distance. Using Rqrb(z = 0) « 
0.5 x lO^Mpc" 3 yr" 1 03l, 

we find that a large number of GRBs, 
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may contribute to the flux at any given time. 

Figure [fl adapted from Rial , presents a comparison of available UHECR data with the 
predictions of a model, where extra-galactic protons in the energy range E p < 10 21 eV 
are produced by cosmologically-distributed sources at a rate and spectrum given by 
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Here, (z) accounts for redshift evolution and <j) (z = 0) = 1. The spectrum above 10 19 eV 
is only weakly dependent on 0(z) since proton energy loss limits their propagation 
distance. For the heavy nuclei component dominating at lower, < 10 19 eV, energy the 
Fly's Eye experimental fit 111711 . dN/dE E~ 3 - 50 , was used. The power-law spectrum 
of accelerated particles, dN/dE oc E~ 2 , has been observed for both non-relativistic 
and relativistic shocks, and is believed to be due to Fermi acceleration in collisionless 
shocks [[H. 

Model predictions are in good agreement with the data of all experiments in the energy 
range 10 19 eV to 10 20 eV (As explained in detail in the various experiments are 
consistent with each other when systematic errors in the absolute energy scale of the 
events are taken into account). The suppression of the flux above ~ 10 19,7 eV is the 
manifestation of the GZK effect. Above 10 20 eV the Fly's Eye, HiRes and Yakutsk 
experiments are in agreement with each other and with the model, while the AGASA 
experiment reports a flux higher by a factor ~ 3. The origin of this discrepancy is unclear. 
The Auger UHECR experiment Il20n . which is currently under construction, is expected 
to significantly reduce the systematic uncertainty in cosmic -ray energy determination 
and to dramatically increase the number of detected UHECRs (thus reducing statistical 
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FIGURE 1. The solid curve shows the energy spectrum derived from the two-component model dis- 
cussed in § 2.2 (with <j)(z) x (1 +z) 3 up to z = 2, following the evolution of star formation rate). The 
dashed curve shows the extra-Galactic component contribution. The "No GZK" curve is an extrapolation 
of the £~ 2 - 75 energy spectrum derived for the energy range of6x 10 18 eV to4x 10 19 eV fl]. Data taken 
from lfl7l[l9ll (AGASA's energy scale was chosen). 



errors). It will provide an accurate determination of the UHECR spectrum to energies 
> 10 20 eV. 

The local UHECR energy production rate, eq. ©, is remarkably similar to the local 
y-ray energy production rate by GRBs, « 10 44 erg Mpc~ 3 yr^ 1 plQ . The similarity of 
the energy production rates was the second main argument suggested in lfl"Dl in support 
of an association between UHECR and GRB sources. 



2.3. GRB model predictions 

The rapid decrease of D c with E p implies, for 5A 1//2 < 10~ 8 G Mpc 1 / 2 and Rqrb < 
1 / Gpc 3 yr, that for some energy E c in the range 10 20 eV < E c < 4x 10 20 eV the number 
of GRBs contributing to the UHECR flux is (on average) 1, Ngrb(> E c ) = 1 (see eq.B}. 
Fig. [2] presents the flux obtained for E c = 1.4 x 10 20 eV in one realization of a Monte- 
Carlo simulation described in [14]. For each realization the positions and times at which 
cosmological GRBs occurred were randomly drawn. Most of the realizations gave an 
overall spectrum similar to that obtained in the realization of Fig. [2] when the brightest 
source of this realization (dominating at 10 20 eV) is not included. At E p < E c , the number 
of sources contributing to the flux is very large, and the UHECR flux received at any 
given time is near the average (the average flux is that obtained when the UHECR 
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FIGURE 2. Results of a Monte-Carlo realization of the bursting sources model, with E c = 1.4 x 
10 20 eV: Thick solid line- overall spectrum in the realization; Thin solid line- average spectrum, this 
curve also gives D c (E p ); Dotted lines- spectra of brightest sources at different energies. 



emissivity is spatially uniform and time independent). At E p > E c , the flux will generally 
be much lower than the average, because there will be no burst within a distance D c (E p ) 
having taken place sufficiently recently. There is, however, a significant probability to 
observe one source with a flux higher than the average. A source similar to the brightest 
one in Fig. [2] appears ~ 5% of the time. 

At any fixed time a given burst is observed in UHECRs only over a narrow range of 
energy, because if a burst is currently observed at some energy E p , then UHECRs of 
much lower energy from this burst have not yet arrived, while higher energy UHECRs 
reached us mostly in the past. Thus, bursting UHECR sources should have narrowly 
peaked energy spectra, and the brightest sources should be different at different energies. 
For steady state sources, on the other hand, the brightest source at high energies should 
also be the brightest one at low energies. The large exposure expected for the Auger 
experiment |20D may allow to distinguish between the two cases, and, if the signature 
of bursting sources is detected, to determine the bursting source properties (e.g. rate, 
energy production). 



3. HIGH ENERGY NEUTRINOS 



3.1. An upper bound 



The energy production rate, eq. [5l sets an upper bound to the neutrino intensity 
produced by sources which, like GRBs and AGN jets, are optically thin for high-energy 
nucleons to py and pp(n) interactions. For sources of this type, the energy generation 
rate of neutrinos can not exceed the energy generation rate implied by assuming that 
all the energy injected as high-energy protons is converted to pions (via py and pp(n) 
interactions). The resulting upper bound (for muon and anti-muon neutrinos, neglecting 
mixing) is |22Q 
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^ z is (a dimensionless parameter) of order unity, which depends on the redshift evolution 
of EpdN p /dE p (see eq. [5]). In order to obtain a conservative upper bound, we adopt 

(EpdN p /dE p ) z= o = 10 44 erg/Mpc 3 yr and a rapid redshift evolution, 4>(z) = (1 +z) 3 up 
to z = 2, following the evolution of star formation rate. This evolution yields ^ ~ 3. The 
WB upper bound is compared in fig. [3] with current experimental limits, and with the 
expected sensitivity of planned neutrino telescopes. The figure indicates that km-scale 
(i.e. giga-ton) neutrino telescopes are needed to detect the expected extra- Galactic flux 
in the energy range of ~ 1 TeV to ~ 1 PeV, and that much larger effective volume is 
required to detect the flux at higher energy. 



3.2. "Generic" 100 TeV fireball neutrinos 

Protons accelerated in the fireball internal shocks, where GRB y-rays are expected 
to be produced, lose energy through photo-production of mesons in interactions with 
fireball photons. The decay of charged pions produced in this interaction results in the 
production of high energy neutrinos. The key relation is between the observed photon 
energy, Ey, and the accelerated proton's energy, E p , at the threshold of the A-resonance. 
In the observer frame, 

£ r £ p = 0.2GeV 2 r 2 . (7) 

For T « 10 2 - 5 and E y = 1 MeV, we see that characteristic proton energies ~ 10 16 eV are 
required to produce pions. Since neutrinos produced by pion decay typically carry 5% 
of the proton energy, production of ~ 10 14 eV neutrinos is expected oh. 

The fraction of energy lost by protons to pions, f K , is f K 0.2 [0,130. Assuming that 
GRBs generate the observed UHECRs, the expected GRB muon and anti-muon neutrino 
flux may be estimated using eq.l6l|pL l22|1 . 

£ 2 4> v « 0.8 x 10" 8 ^GeVcm _2 s _1 sr _1 . (8) 
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FIGURE 3. The upper bound on the extra-Galactic muon and tau neutrino flux (lower-curve: no source 
evolution, upper curve: source evolution following star formation rate), assuming oscillations lead to 
<P Ve : 4> VfI : < I ) Vt =1:1 : 1, compared with experimental upper bounds (solid lines) of optical Cerenkov 
experiments (BAIKAL Q, AMANDA Q), and of coherent Cerenkov radio experiments (RICE lEll . 
GLUE 12^1 ; see also ll27ll ). The curve labelled "GZK" shows the intensity due to interaction with 
micro-wave background photons. Dashed curves show the expected sensitivity of 0.1 Gton (AMANDA, 
ANTARES, NESTOR) and 1 Gton (IceCube, NEMO) optical Cerenkov detectors [28], of the coherent 
radio Cerenkov (balloon) experiment ANITA [27] and of the Auger air-shower detector (sensitivity to v T ) 
lEsill . Space air-shower detectors (OWL-AIRWATCH) may also achieve the sensitivity required to detect 
fluxes lower than the WB bound at energies > 10 18 eV 112911 . 



This neutrino spectrum extends to ~ 10 eV, and is suppressed at higher energy due 
to energy loss of pions and muons flHIli. Eq.i implies a detection rate of ~ 20 
neutrino-induced muon events per year (over An sr) in a cubic -km detector 

flUD]. Since 

GRB neutrino events are correlated both in time and in direction with gamma-rays, their 
detection is practically background free. 



3.3. TeV neutrinos 

The 100 TeV neutrinos discussed in the previous section are produced at the same 
region where GRB y-rays are produced. Their production is a generic prediction of 
the fireball model. It is a direct consequence of the assumptions that energy is carried 
from the underlying engine as kinetic energy of protons and that y-rays are produced by 
synchrotron emission of shock accelerated particles. Neutrinos may be produced also in 
other stages of fireball evolution, at energies different than 100 TeV. The production of 
these neutrinos is dependent on additional model assumptions. We discuss below some 
examples related to the GRB progenitor. For a more detailed discussion see [32, 1] and 
references therein. 

The most widely discussed progenitor scenarios for long-duration GRBs involve core 
collapse of massive stars. In these "collapsar" models, a relativistic jet breaks through the 
stellar envelope to produce a GRB. For extended or slowly rotating stars, the jet may be 
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FIGURE 4. Flavor and anti-particle content of the flux of astrophysical neutrinos produced by pion 
decay, for pion energy spectrum (at production) dn jdE % °= E^ 2 and electromagnetic pion energy loss 
rate dE K /dt °= E\, as expected for internal shocks in GRB fireballs. Left: The ratio between ^v^fvr) 
and <t>v e (solid line), with 90% CL lines of (dashed) and v T (dotted) fluxes (here, <J> V/ stands for the 
combined flux of V/ and V/). Right: The ratio of v e to total v flux on Earth, solid (dashed) line for neutrinos 
produced by py (pp) interactions. Eq „ is the muon energy for which the muon life time is comparable to 
its electromagnetic energy loss time. Eq.^ ~ 10 3 TeV for internal shocks in GRB fireballs. 



unable to break through the envelope. Both penetrating (GRB producing) and "choked" 
jets can produce a burst of ~ 5 TeV neutrinos by interaction of accelerated protons with 
jet photons, while the jet propagates in the envelope 11331. 13411 . The estimated event rates 
may exceed ~ 10 2 events per yr in a km- scale detector, depending on the ratio of non- 
visible to visible fireballs. A clear detection of non-visible GRBs with neutrinos may be 
difficult due to the low energy resolution for muon-neutrino events, unless the associated 
supernova photons are detected. In the two-step "supranova" model, interaction of the 
GRB blast wave with the supernova shell can lead to detectable neutrino emission, either 
through nuclear collisions with the dense supernova shell or through interaction with the 
intense supernova and backscattered radiation field ([351] . 



3.4. Energy dependent neutrino flavor ratio 

Although a 1 : 1 : 1 flavor ratio appears to be a robust prediction of models where 
neutrinos are produced by pion decay, energy dependence of the flavor ratio is a generic 
feature of models of high energy astrophysical neutrino sources Q36Q. Pions are typically 
produced in environments where they may suffer significant energy losses prior to decay, 
due to interaction with radiation and magnetic fields [0,[3(1]. Since the pion life time is 
shorter than the muon's, at sufficiently high energy the probability for pion decay prior 
to significant energy loss is higher than the corresponding probability for muon decay. 
This leads to suppression at high energy of the relative contribution of muon decay to the 
neutrino flux. The flavor ratio is modified to : 1 : at the source, implying 1 : 1.8 : 1.8 
ratio on Earth 113611 . 

Figure |4] presents the expected energy dependence of flavor and anti-particle content 



for internal shocks in GRB fireballs. Since the transition is expected at ~ 100 TeV, it 
may be detected by km- scale v telescopes. 
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